Stochastic models for regulatory networks of the genetic toggle switch

Bistability arises within a wide range of biological systems from the A phage switch in bacteria to cellular signal transduction pathways in mammalian cells. Changes in regulatory mechanisms may result in genetic switching in a bistable system. Recently, more and more experimental evidence in the form of bimodal population distributions indicates that noise plays a very important role in the switching of bistable systems. Although deterministic models have been used for studying the existence of bistability properties under various system conditions, these models cannot realize cell-to-cell fluctuations in genetic switching. However, there is a lag in the development of stochastic models for studying the impact of noise in bistable systems because of the lack of detailed knowledge of biochemical reactions, kinetic rates, and molecular numbers. in this work, we develop a previously undescribed general technique for developing quantitative stochastic models for large-scale genetic regulatory networks by introducing Poisson random variables into deterministic models described by ordinary differential equations. Two stochastic models have been proposed for the genetic toggle switch interfaced with either the SOS signaling pathway or a quorum-sensing signaling pathway, and we have successfully realized experimental results showing bimodal population distributions. Because the introduced stochastic models are based on widely used ordinary differential equation models, the success of this work suggests that this approach is a very promising one for studying noise in large-scale genetic regulatory networks.

A genetic two-factor model of the covariation among a subset of Multidimensional Aptitude Battery and Wechsler Adult Intelligence Scale - Revised subtests

The phenotypic and genetic factor structure of performance on five Multidimensional Aptitude Battery (MAB) subtests and one Wechsler Adult Intelligence Scale-Revised (WAIS-R) subtest was explored in 390 adolescent twin pairs (184 monozygotic [MZ]; 206 dizygotic (DZ)). The temporal stability of these measures was derived from a subsample of 49 twin pairs, with test-retest correlations ranging from .67 to .85. A phenotypic factor model, in which performance and verbal factors were correlated, provided a good fit to the data. Genetic modeling was based on the phenotypic factor structure, but also took into account the additive genetic (A), common environmental (C), and unique environmental (E) parameters derived from a fully saturated ACE model. The best fitting model was characterized by a genetic correlated two-factor structure with specific effects, a general common environmental factor, and overlapping unique environmental effects. Results are compared to multivariate genetic models reported in children and adults, with the most notable difference being the growing importance of common genes influencing diverse abilities in adolescence. (C) 2003 Elsevier Inc. All rights reserved.

Using identity-by-de scent information in affected sib pairs to increase the efficiency of genetic association studies

T he aim of this study was to determine whether identity-by-descent (IBD) information for affected sib pairs (ASPs) can be used to select a sample of cases for a genetic case-control study which will provide more power for detecting association with loci in a known linkage region. By modeling the expected frequency of the disease allele in ASPs showing IBD sharing of 0, 1, or 2 alleles, and considering additive, recessive, and dominant disease models, we show that cases selected from IBD 2 families are best for this purpose, followed by those selected from IBD 1 families; least useful are cases selected from IBD 0 families.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Do the Genetic or Environmental Determinants of Anxiety and Depression Change with Age? A Longitudinal Study of Australian Twins

Because the determinants of anxiety and depression in late adolescence and early adulthood may differ from those in later life, we investigated the temporal stability and magnitude of genetic and environmental correlates of symptoms of anxiety and depression across the life span. Data were collected from a population-based Australian sample of 4364 complete twin pairs and 777 singletons aged 20 to 96 years who were followed-up over three studies between 1980 and 1996. Each study contained the 14-item self-report DSSI/sAD scale which was used to measure recently experienced symptoms of anxiety and depression. Symptom scores were then divided and assigned to age intervals according to each subject's age at time of participation. We fitted genetic simplex models to take into account the longitudinal nature of the data. For male anxiety and depression, the best fitting simplex models comprised a single genetic innovation at age 20 which was transmitted, and explained genetic variation in anxiety and depression at ages 30, 40, 50 and 60. Most of the lifetime genetic variation in female anxiety and depression could also be explained by innovations at age 20 which were transmitted to all other ages; however, there were also smaller age-dependent genetic innovations at 30 for anxiety and at 40 and 70 for depression. Although the genetic determinants of anxiety and depression appear relatively stable across the life-span for males and females, there is some evidence to support additional mid-life and late age gene action in females for depression. The fact that mid-life onset for anxiety occurs one decade before depression is also consistent with a causal relationship (anxiety leading to depression) between these conditions. These findings have significance for large scale depression prevention projects.

Genetic Epidemiology of Alcohol-Induced Blackouts

Background: Alcohol-induced blackouts (ie, periods of anterograde amnesia) have received limited recent research attention. Objective: To examine the genetic epidemiology of lifetime blackouts and having had 3 or more blackouts in a year, including analyses controlling for the frequency of intoxication. Design, Setting, and Participants: Members of the young adult Australian Twin Register, a volunteer twin panel born between January 1, 1964, and December 3 1, 1971, were initially registered with the panel as children by their parents between 1980 and 1982. They underwent structured psychiatric telephone inter-views from February 1996 through September 2000. The current sample contains 2324 monozygotic and dizygotic twin pairs (mean [SDI age 29.9 [2.5] years) for whom both twins' responses were coded for blackout questions and for frequency of intoxication. Main Outcome Measure: Data on lifetime blackouts and having had 3 or more blackouts in a year were collected within an examination of the genetic epidemiology of alcoholism. Results: A lifetime history of blackouts was reported by 39.3% of women and 52.4% of men; 11.4% of women and 20.9% of men reported having had 3 or more blackouts in a year. The heritability of lifetime blackouts was 52.5% and that of having had 3 or more blackouts in a year was 57.8%. Models that controlled for frequency of intoxication found evidence of substantial genetic contribution unique to risk for the blackouts and a significant component of genetic risk shared with frequency of intoxication. Conclusions: The finding of a substantial genetic contribution to liability for alcohol-induced blackouts including a component of genetic loading shared with frequency of intoxication may offer important additional avenues to investigate susceptibility to alcohol-related problems.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Genetic effects on alcohol dependence risk: re-evaluating the importance of psychiatric and other heritable risk factors

Background. Genetic influences have been shown to play a major role in determining the risk of alcohol dependence (AD) in both women and men; however, little attention has been directed to identifying the major sources of genetic variation in AD risk. Method. Diagnostic telephone interview data from young adult Australian twin pairs born between 1964 and 1971 were analyzed. Cox regression models were fitted to interview data from a total of 2708 complete twin pairs (690 MZ female, 485 MZ male, 500 DZ female, 384 DZ male, and 649 DZ female/male pairs). Structural equation models were fitted to determine the extent of residual genetic and environmental influences on AD risk while controlling for effects of sociodemographic and psychiatric predictors on risk. Results. Risk of AD was increased in males, in Roman Catholics, in those reporting a history of major depression, social anxiety problems, and conduct disorder, or (in females only) a history of suicide attempt and childhood sexual abuse; but was decreased in those reporting Baptist, Methodist, or Orthodox religion, in those who reported weekly church attendance, and in university-educated males. After allowing for the effects of sociodemographic and psychiatric predictors, 47 % (95 % CI 28-55) of the residual variance in alcoholism risk was attributable to additive genetic effects, 0 % (95 % CI 0-14) to shared environmental factors, and 53 % (95 % CI 45-63) to non-shared environmental influences. Conclusions. Controlling for other risk factors, substantial residual heritability of AD was observed, suggesting that psychiatric and other risk factors play a minor role in the inheritance of AD.

Multivariate quantitative genetics and the Lek Paradox: Genetic variance in male sexually selected traits of Drosophila serrata under field control

Single male sexually selected traits have been found to exhibit substantial genetic variance, even though natural and sexual selection are predicted to deplete genetic variance in these traits. We tested whether genetic variance in multiple male display traits of Drosophila serrata was maintained under field conditions. A breeding design involving 300 field-reared males and their laboratory-reared offspring allowed the estimation of the genetic variance-covariance matrix for six male cuticular hydrocarbons (CHCs) under field conditions. Despite individual CHCs displaying substantial genetic variance under field conditions, the vast majority of genetic variance in CHCs was not closely associated with the direction of sexual selection measured on field phenotypes. Relative concentrations of three CHCs correlated positively with body size in the field, but not under laboratory conditions, suggesting condition-dependent expression of CHCs under field conditions. Therefore condition dependence may not maintain genetic variance in preferred combinations of male CHCs under field conditions, suggesting that the large mutational target supplied by the evolution of condition dependence may not provide a solution to the lek paradox in this species. Sustained sexual selection may be adequate to deplete genetic variance in the direction of selection, perhaps as a consequence of the low rate of favorable mutations expected in multiple trait systems.

On dendrites in Down syndrome and DS murine models: a spiny way to learn

Since the discovery in the 1970s that dendritic abnormalities in cortical pyramidal neurons are the most consistent pathologic correlate of mental retardation, research has focused on how dendritic alterations are related to reduced intellectual ability. Due in part to obvious ethical problems and in part to the lack of fruitful methods to study neuronal circuitry in the human cortex, there is little data about the microanatomical contribution to mental retardation. The recent identification of the genetic bases of some mental retardation associated alterations, coupled with the technology to create transgenic animal models and the introduction of powerful sophisticated tools in the field of microanatomy, has led to a growth in the studies of the alterations of pyramidal cell morphology in these disorders. Studies of individuals with Down syndrome, the most frequent genetic disorder leading to mental retardation, allow the analysis of the relationships between cognition, genotype and brain microanatomy. In Down syndrome the crucial question is to define the mechanisms by which an excess of normal gene products, in interaction with the environment, directs and constrains neural maturation, and how this abnormal development translates into cognition and behaviour. In the present article we discuss mainly Down syndrome-associated dendritic abnormalities and plasticity and the role of animal models in these studies. We believe that through the further development of such approaches, the study of the microanatomical substrates of mental retardation will contribute significantly to our understanding of the mechanisms underlying human brain disorders associated with mental retardation. (C) 2004 Elsevier Ltd. All rights reserved.

Bistability and switching in the lysis/lysogeny genetic regulatory network of bacteriophage lambda

Bistability and switching are two important aspects of the genetic regulatory network of phage. Positive and negative feedbacks are key regulatory mechanisms in this network. By the introduction of threshold values, the developmental pathway of A phage is divided into different stages. If the protein level reaches a threshold value, positive or negative feedback will be effective and regulate the process of development. Using this regulatory mechanism, we present a quantitative model to realize bistability and switching of phage based on experimental data. This model gives descriptions of decisive mechanisms for different pathways in induction. A stochastic model is also introduced for describing statistical properties of switching in induction. A stochastic degradation rate is used to represent intrinsic noise in induction for switching the system from the lysogenic pathway to the lysis pathway. The approach in this paper represents an attempt to describe the regulatory mechanism in genetic regulatory network under the influence of intrinsic noise in the framework of continuous models. (C) 2003 Elsevier Ltd. All rights reserved.

A genetic estimation algorithm for parameters of stochastic ordinary differential equations

A generic method for the estimation of parameters for Stochastic Ordinary Differential Equations (SODEs) is introduced and developed. This algorithm, called the GePERs method, utilises a genetic optimisation algorithm to minimise a stochastic objective function based on the Kolmogorov-Smirnov statistic. Numerical simulations are utilised to form the KS statistic. Further, the examination of some of the factors that improve the precision of the estimates is conducted. This method is used to estimate parameters of diffusion equations and jump-diffusion equations. It is also applied to the problem of model selection for the Queensland electricity market. (C) 2003 Elsevier B.V. All rights reserved.

Genetic time-series analysis identifies a major QTL for in vivo alcohol metabolism not predicted by in vitro studies of structural protein polymorphism at the ADH1B or ADH1C loci

After ingestion of a standardized dose of ethanol, alcohol concentrations were assessed, over 3.5 hours from blood (six readings) and breath (10 readings) in a sample of 412 MZ and DZ twins who took part in an Alcohol Challenge Twin Study (ACTS). Nearly all participants were subsequently genotyped on two polymorphic SNPs in the ADH1B and ADH1C loci known to affect in vitro ADH activity. In the DZ pairs, 14 microsatellite markers covering a 20.5 cM region on chromosome 4 that includes the ADH gene family were assessed, Variation in the timed series of autocorrelated blood and breath alcohol readings was studied using a bivariate simplex design. The contribution of a quantitative trait locus (QTL) or QTL's linked to the ADH region was estimated via a mixture of likelihoods weighted by identity-by-descent probabilities. The effects of allelic substitution at the ADH1B and ADH1C loci were estimated in the means part of the model simultaneously with the effects sex and age. There was a major contribution to variance in alcohol metabolism due to a QTL which accounted for about 64% of the additive genetic covariation common to both blood and breath alcohol readings at the first time point. No effects of the ADH1B*47His or ADH1C*349Ile alleles on in vivo metabolism were observed, although these have been shown to have major effects in vitro. This implies that there is a major determinant of variation for in vivo alcohol metabolism in the ADH region that is not accounted for by these polymorphisms. Earlier analyses of these data suggested that alcohol metabolism is related to drinking behavior and imply that this QTL may be protective against alcohol dependence.